Cosmography of f(R) gravity 
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It is nowadays accepted that the universe is undergoing a phase of accelerated expansion as tested 
by the Hubble diagram of Type la Supernovae (SNela) and several LSS observations. Future SNela 
surveys and other probes will make it possible to better characterize the dynamical state of the 
universe renewing the interest in cosmography which allows a model independent analysis of the 
distance - redshift relation. On the other hand, fourth order theories of gravity, also referred to as 
f(R) gravity, have attracted a lot of interest since they could be able to explain the accelerated 
expansion without any dark energy. We show here how it is possible to relate the cosmographic 
parameters (namely the deceleration qo, the jerk jo, the snap so and the lerk lo parameters) to the 
present day values of f(R) and its derivatives f^ n '(R) — d n f /dR n (with n = 1, 2, 3) thus offering a 
new tool to constrain such higher order models. Our analysis thus offers the possibility to relate the 
model independent results coming from cosmography to the theoretically motivated assumptions of 
f(R) cosmology. 
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I. INTRODUCTION 



As soon as astrophysicists realized that Type la Su- 
pernovae (SNela) were standard candles, it appeared ev- 
ident that their high luminosity should make it possi- 
ble to build a Hubble diagram, i.e. a plot of the dis- 
tance - redshift relation, over some cosmologically inter- 
esting distance ranges. Motivated by this attractive con- 
sideration, two independent teams started SNela surveys 
leading to the unexpected discovery that the universe ex- 
pansion is speeding up rather than decelerating [l| . This 
surprising result has now been strengthened by more re- 
cent data coming from SNela surveys 0, H, 0, [H, @, 0] , 
large scale structure Q and cosmic microwave back- 
ground (CMBR) anisotropy spectrum @, G3 G3 . This 
large dataset coherently points toward the picture of a 
spatially flat universe undergoing an accelerated expan- 
sion driven by a dominant negative pressure fluid, typi- 
cally referred to as dark energy (l2| . 

While there is a wide consensus on the above scenario 
depicted by such good quality data, there is a similarly 
wide range of contrasting proposals to solve the dark 
energy puzzle. Surprisingly, the simplest explanation, 
namely the cosmological constant A Jl3f , is also the best 
one from a statistical point of view [l4j . Unfortunately, 
the well known coincidence and 120 orders of magnitude 
problems render A a rather unattractive solution from a 
theoretical point of view. Inspired by the analogy with 
inflation, a scalar field (f>, dubbed quintessence [15] , has 
then been proposed to give a dynamical A term in or- 
der to both fit the data and avoid the above problems. 
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However, such models are still plagued by difficulties on 
their own, such as the almost complete freedom in the 
choice of the scalar field potential and the fine tuning of 
the initial conditions. Needless to say, a plethora of al- 
ternative models are now on the market all sharing the 
main property to be in agreement with observations, but 
relying on completely different physics. 

Notwithstanding their differences, all the dark energy 
based theories assume that the observed acceleration is 
the outcome of the action of an up to now undetected 
ingredient to be added to the cosmic pie. In terms of the 
Einstein equations, G M „ = xF^v, such models are simply 
modifying the right hand side including in the stress - 
energy tensor something more than the usual matter and 
radiation components. 

As a radically different approach, one can also try to 
leave unchanged the source side, but rather modifying 
the left hand side. In a sense, one is therefore interpret- 
ing cosmic speed up as a first signal of the breakdown of 
the laws of physics as described by the standard General 
Relativity (GR). Since this theory has been experimen- 
tally tested only up to the Solar System scale, there is 
no a priori theoretical motivation to extend its validity 
to extraordinarily larger scales such as the cosmologi- 
cal ones (e.g. the last scattering surface!). Extending 
GR, not giving up to its positive results, opens the way 
to a large class of alternative theories of gravity rang- 
ing from extra - dimensions [161 ] to nonminimally coupled 
scalar fields [13, EH- In particular, we will be interested 
here in fourth order theories [ijj H(| based on replac- 
ing the scalar curvature R in the Hilbert-Einstein action 
with a generic analytic function f{R) which should be re- 
constructed starting from data and physically motivated 
issues. Also referred to as f{R) gravity, these models 
have been shown to be able to both fit the cosmological 
data and evade the Solar System constraints in several 
physically interesting cases (Hi, [H, SI [H, Hi . 
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It is worth noting that both dark energy models and 
modified gravity theories have shown to be in agreement 
with the data. As a consequence, unless higher precision 
probes of the expansion rate and the growth of structure 
will be available, these two rival approaches could not 
be discriminated. This confusion about the theoretical 
background suggests that a more conservative approach 
to the problem of cosmic acceleration, relying on as less 
model dependent quantities as possible, is welcome. A 
possible solution could be to come back to the cosmog- 
raphy 26] rather than finding out solutions of the Fried- 
mann equations and testing them. Being only related to 
the derivatives of the scale factor, the cosmographic pa- 
rameters make it possible to fit the data on the distance - 
redshift relation without any a priori assumption on the 
underlying cosmological model: in this case, the only as- 
sumption is that the metric is the Robertson - Walker one 
(and hence not relying on the solution of cosmic equa- 
tions). Almost a century after Hubble discovery of the 
expansion of the universe, we could now extend cosmog- 
raphy beyond the search for the value of the Hubble con- 
stant. The SNela Hubble diagram extends up to z = 1.7 
thus invoking the need for, at least, a fifth order Taylor 
expansion of the scale factor in order to give a reliable 
approximation of the distance - redshift relation. As a 
consequence, it could be, in principle, possible to esti- 
mate up to five cosmographic parameters, although the 
still too small dataset available does not allow to get a 
precise and realistic determination of all of them. 

Once these quantities have been determined, one could 
use them to put constraints on the models. In a sense, 
we are reverting the usual approach consisting in deriving 
the cosmographic parameters as a sort of byproduct of an 
assumed theory. Here, we follow the other way around 
expressing the model characterizing quantities as a func- 
tion of the cosmographic parameters. Such a program is 
particularly suited for the study of fourth order theories 
of gravity. As is well known, the mathematical difficul- 
ties entering the solution of fourth order field equations 
make it quite problematic to find out analytical expres- 
sions for the scale factor and hence predict the values of 
the cosmographic parameters. A key role in /(i?) grav- 
ity is played by the choice of the f{R) function. Under 
quite general hypotheses, we will derive useful relations 
among the cosmographic parameters and the present day 
value of /(") (R) = d n f/dR n , with n = 0, . . . , 3, whatever 
f(R) is 1 . Once the cosmographic parameters will be de- 
termined, this method will allow us to investigate the 
cosmography of f(R) theories. 

The layout of the paper is as follows. Sects. II and III 
are devoted to introducing the basic notions of the cosmo- 
graphic parameters and f(R) gravity respectively, sum- 



1 As an important remark, we stress that our derivation will rely 
on the metric formulation of f(R) theories, while we refer the 
reader to |27|| for a similar work in the Palatini approach. 



marizing the main formulae we will use later. Sect. IV 
contains the main result of the paper demonstrating how 
the f(R) derivatives can be related to the cosmographic 
parameters. Since these latter are not well determined 
today, we will discuss, in Sect. V, how these formulae 
can be adapted to a different parameterization relying 
on expressing the cosmographic parameters in terms of a 
phenomenological assumption for the dark energy equa- 
tion of state. Sect. VI illustrates a possible application 
of the relation among f(R) derivatives and cosmographic 
parameters showing how one can constrain the parame- 
ters of a given f(R) model. Since future data will likely 
determine with a sufficient precision at least the first two 
cosmographic parameters, it is worth estimating how this 
will impact on the determination of the f(R) quantities, 
which is the argument of Sect. VII. We then summarize 
and conclude in Sect. VIII. 



II. COSMOGRAPHIC PARAMETERS 

Standard candles (such as SNela and, to a limited ex- 
tent, gamma ray bursts) are ideal tools in modern cos- 
mology since they make it possible to reconstruct the 
Hubble diagram, i.e. the redshift - distance relation up 
to high redshift values. It is then customary to assume 
a parameterized model (such as the concordance ACDM 
one, or any other kind of dark energy scenario) and con- 
trasting it against the data to check its viability and con- 
straints its characterizing parameters. As it is clear, such 
an approach is model dependent so that some doubts al- 
ways remain on the validity of the constraints on derived 
quantities as the present day values of the deceleration 
parameter and the age of the universe. In order to over- 
come such a problem, one may resort to cosmography, 
i.e. expanding the scale factor in Taylor series with re- 
spect to the cosmic time [26j . Such an expansions leads 
to a distance - redshit relation which only relies on the 
assumption of the Robertson - Walker metric thus being 
fully model independent since it does not depend on the 
particular form of the solution of cosmic equations. To 
this aim, it is convenient to introduce the following func- 
tions HE [H : 



a at 

_ 1 d 2 a , 

q ~ "TTT H 

a dt z 

a dt A 
a dt* 
a dt 5 
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which are usually referred to as the Hubble, deceleration, 
jerk, snap and lerk parameters [29|, respectively 2 . Their 
present day values (which we will denote with a subscript 
0) may be used to characterize the evolutionary status of 
the Universe. For instance, qo < denotes an acceler- 
ated expansion, while jo allows to discriminate among 
different accelerating models. 

It is then a matter of algebra to demonstrate the fol- 
lowing useful relations : 



2% up to z ~ 2 over a wide range of the CPL parameter 
space. On the other hand, halting the Taylor expansion 
to a lower order may introduce significant deviation for 
z > 1 that can potentially bias the analysis if the mea- 
surement errors are as small as those predicted for future 
SNela surveys. We are therefore confident that our fifth 
order expansion is both sufficient to get an accurate dis- 
tance modulus over the redshift range probed by SNela 
and necessary to avoid dangerous biases. 



H=-H 2 (l + q) 



(2) 



III. f(R) GRAVITY 



H = H 3 (j +3q + 2) 



H = H i [s - 4j - 3q(q + 4) - 6] 



(3) 



(4) 



Much interest has been recently devoted to a form of 
quintessence induced by curvature according to which the 
present universe is filled by pressureless dust matter only 
and the acceleration is the result of the modified Fried- 
mann equations obtained by replacing the Ricci curva- 
ture scalar R with a generic function /(i?) in the gravity 
action (l9l . [20| . Under the assumption of a flat universe, 
the Hubble parameter is therefore determined by 4 : 



d 4 H/dt A = H 5 [l- 5s + 10(q + 2)j + 30(q + 2)q + 24] , 

(5) 

where a dot denotes derivative with respect to the cosmic 
time t. Eqs.@ - make it possible to relate the deriva- 
tive of the Hubble parameter to the other cosmographic 
parameters. The distance - redshift relation may then be 
obtained starting from the Taylor expansion of a(t) along 
the lines described in [2^, [3l|, [12] • The result for the fifth 
order is reported in Appendix A. 

It is worth stressing that the definition of the cosmo- 
graphic parameters only relies on the assumption of the 
Robertson - Walker metric. As such, it is however difficult 
to state a priori to what extent the fifth order expansion 
provides an accurate enough description of the quanti- 
ties of interest. Actually, the number of cosmographic 
parameters to be used depends on the problem one is in- 
terested in. As we will see later, we are here concerned 
only with the SNela Hubble diagram so that we have to 
check that the distance modulus /i cp (z) obtained using 
the fifth order expansion of the scale factor is the same 
(within the errors) as the one hde(z) of the underlying 
physical model. Being such a model of course unknown, 
one can adopt a phenomenological parameterization for 
the dark energy 3 EoS and look at the percentage devi- 
ation A/i//i£>£ as function of the EoS parameters. We 
have carried out such exercise using the CPL model in- 
troduced later and verified that A/i/ /ide is an increasing 
function of z (as expected) , but still remains smaller than 



f'(R) 



(G) 



where the prime denotes derivative with respect to R 
and Pcurv is the energy density of an effective curvature 
fluid 5 : 



lf(R) - Rf'(R)] - 3HRf"(R) 



f'(R) 



(7) 

Assuming there is no interaction between the matter 
and the curvature terms (we are in the so-called Jordan 
frame), the matter continuity equation gives the usual 
scaling pM = pAi(t = t Q )a~ 3 = 3HqHmcl~ 3 , with Am 
the present day matter density parameter. The continu- 
ity equation for p curv then reads : 



3H(l 



>)Pc 



3H^n M Rf"(R) _ 3 , , 

9 a \°) 

lf'(R)} 



with 



= -1 



Rf"(R) + R Rf"'{R) - Hf"(R) 



{f(R) - Rf'(R)] /2 - ZHRf"{R) 



(9) 



the barotropic factor of the curvature fluid. It is worth 
noticing that the curvature fluid quantities p CU rv and 



2 Note that the use of the jerk parameter to discriminate between 
different models was also proposed in |3d| in the context of the 
statefinder parametrization. 

3 Note that one can always use a phenomenological dark energy 
model to get a reliable estimate of the scale factor evolution even 
if the correct model is a fourth order one. 



4 We use here natural units such that 8-7rG = 1. 

5 Note that the name curvature fluid does not refer to the FRW 
curvature parameter k, but only takes into account that such a 
term is a geometrical one related to the scalar curvature R. 
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w curv only depends on f(R) and its derivatives up to 
the third order. As a consequence, considering only their 
present day values (which may be naively obtained by re- 
placing R with i?o everywhere), two f(R) theories shar- 
ing the same values of f(Ro), f'(Ro), f"{R ), f"(R ) 
will be degenerate from this point of view 6 . 

Combining Eq.© with Eq.©, one finally gets the fol- 
lowing master equation for the Hubble parameter : 



H = - 



2f'(R) 



+R 



Rf"'(R) - Hf" (#)] } 



■Rf"(R)- 



(10) 



Expressing the scalar curvature R as function of the Hub- 
ble parameter as : 



R 



6[H + 2H' 



(11) 



and inserting the result into Eq. (fTU)) . one ends with 
a fourth order nonlinear differential equation for the 
scale factor a(t) that cannot be easily solved also for 
the simplest cases (for instance, f(R) oc R n ). More- 
over, although technically feasible, a numerical solution 
of Eq.(flQ|) is plagued by the large uncertainties on the 
boundary conditions (i.e., the present day values of the 
scale factor and its derivatives up to the third order) that 
have to be set to find out the scale factor. 



IV. f(R) DERIVATIVES VS COSMOGRAPHY 

Motivated by these difficulties, we approach now the 
problem from a different viewpoint. Rather than choos- 
ing a parameterized expression for /(i?) and then numer- 
ically solving Eq.(Un]) for given values of the boundary 
conditions, we try to relate the present day values of its 
derivatives to the cosmographic parameters (qo,jo, sq,Iq) 
so that constraining them in a model independent way 
gives us a hint for what kind of f(R) theory could be 
able to fit the observed Hubble diagram 7 . 

As a preliminary step, it is worth considering again the 
constraint equation (jlip . Differentiating with respect to 
t, we easily get the following relations : 



6 One can argue that this is not strictly true since different f(R) 
theories will lead to different expansion rate H(t) and hence dif- 
ferent present day values of R and its derivatives. However, it 
is likely that two f(R) functions that exactly match each other 
up to the third order derivative today will give rise to the same 
H(t) at least for t ~ to so that (Rq,Ro,Ro) will be almost the 
same. 

7 Note that a similar analysis, but in the context of the energy 
conditions in f(R), has yet been presented in pjjjj ]. However, 
in that paper, the author give an expression for f(R) and then 
compute the snap parameter to be compared to the observed one. 
On the contrary, our analysis does not depend on any assumed 
functional expression for f(R). 



R = -6 (^H + AHH 
R = -&(h + AHH + 4ff 2 



(12) 



R = -6 



(d^H/dt 4 



AHH + 12HH 



Evaluating these at the present time and using Eqs.(TJ|)- 
([5]) , one finally gets : 



Ro = -6ff 2 (l - q Q ) , 



Ro = -6if 3 (jo - q Q - 2) , 



Ro = -6H 4 (s + q 2 a 



6) , 



(13) 



(14) 



(15) 



R 



-6ff 5 [J„ -s + 2{q + 4)jo - 6(3 9o + 8)q 



24] , 
(16) 

which will turn out to be useful in the following. 

Let us now come back to the expansion rate and master 
equations (|6|) and (|10|). Since they have to hold along 
the full evolutionary history of the universe, they naively 
hold also at the present day. As a consequence, we may 
evaluate them in t = to thus easily obtaining : 



m = 



f'(Ro) 

f(Ro) - Rof'(R ) 



6H R f"(R 



6/'(i?o) 



(17) 



H 



1 »Af 



2/W 

+ (Ro 



HqRo ) f"(Ro) 



V'(Ro) 



(18) 



Using Eqs.(J2]) - (JI) and (fT3| - (fl6|) . we can rearrange 
Eqs. (fi"7|) and (fT8|) as two relations among the Hubble con- 
stant Hq and the cosmographic parameters (qo, jo, So), 
on one hand, and the present day values of f(R) and 
its derivatives up to third order. However, two further 
relations are needed in order to close the system and 
determine the four unknown quantities f(Ro), f'(Ro), 
f"(Ro), f"'(Ro). A first one may be easily obtained by 
noting that, inserting back the physical units, the rate 
expansion equation reads : 



H 2 



8nG 
3f'(R) 



+ P curv 
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which clearly shows that, in f(R) gravity, the Newtonian 
gravitational constant G is replaced by an effective (time 
dependent) G e // = G/f'(R). On the other hand, it is 
reasonable to assume that the present day value of G e / / 
is the same as the Newtonian one so that we get the 
simple constraint : 



G eff (z = 0) = G^f'(Ro) = l 



(19) 



In order to get the fourth relation we need to close the 
system, we first differentiate both sides of Eq. pTj|) with 
respect to t. We thus get : 



H 



R 2 f"'(R) + (R- HRj f"(R) + 3H$n M a- 



-l 



Rf"(R)\ [f'(R)Y 
R?f(™)(R) + (zRR - HR 2 ^ f"'(R) 

mR) 

Ck' - HR + HRj f"(R) - 9H%n M Ha- 
2f'(R) 



,(20) 



with fW(R) = d i f/dR i . Let us now suppose that f(R) 
may be well approximated by its third order Taylor ex- 
pansion in R — Rq, i.e. we set : 



f(R) = f(Ro) + f , (R a )(R-Ro) + lf"(R )(R-Ro) 2 



C» 



f'"(R )(R-Ro) 3 



(21) 



In such an approximation, it is f( n '(R) = d n f/R n = 
for n > 4 so that naively f^(Ro) = 0. Evaluating then 
Eq. ([2H)l at the present day, we get : 



H 



Rir(Ro) + [Ro- H Q R Q f"(R ) + 3H 2 n 



I'M 



2 [R f"(Ro)\ [f'(Ro)] 2 
^3R Rq - HRlj f'"(R ) 

WW) 

(^Ro — HqRo + HqRq^J f"(R ) — 9i?o^M 
2f'(Ro) 



(22) 



We can now schematically proceed as follows. Evaluate 
Eqs.([5]) - (O at z — and plug these relations into the 
left hand sides of Eqs.(TF|), dHJ), ([22]). Insert Eqs.([13])- 
(TiTj]) into the right hand sides of these same equations 
so that only the cosmographic parameters (<?o> jo>SO)fo) 
and the f(R) related quantities enter both sides of these 
relations. Finally, solve them under the constraint (fl"9|) 
with respect to the present day values of f{R) and its 
derivatives up to the third order. After some algebra, 
one ends up with the desired result : 



/Cgo) 
6H 2 



~Pa(qn,jo,s a ,lo)tlM + Qo(go, jo, s ,l ) 



f"(Ro) 



(23) 

f(Ro) = 1 , (24) 

TMgo, jo,so)^M + Q 2 (go, jo, so) 
fc(qoJo,so,k) 

"Pa(qo,jo,so,lo)£lM + Q.3(qo,jo,so,k) 



(25) 



(jo - ?o - 2)U(q , j , s , h) 



where we have defined : 



(26) 



Vo = (jo - qo - 2)l 

- (3s + 7 j + 6q 2 + 41q + 22)s 

- j(3g + 16)jo + 20q 2 + 64q + 12] j 

- (3^ + 25q 3 + 96q 2 + 72q + 20) , (27) 

Qo = (<7o - jo9o + 2q )l 

+ [3q s + {Aqo + 6)j + 6g§ + Uq 2 + 22q - 12] s 

+ [2j 2 + (3q 2 + 10q - 6)j + 17q$ + 52q 2 + 54q 

+ 36] j + 3g§ + 28^ + 1 18q 3 + 72q 2 - 76q 

- 64 , (28) 



V 2 = 9s + 6j + 9q 2 + 66q + 42 

Q 2 = -{6(<7o + l)s 
+ [2j - 2(1 - (7o)bo 
+ 6^ + 50^ + 74g + 32} . 



(29) 



(30) 



V 3 = 3l + 3s - % + 4)j - (45gg + 78g + 12) , (31) 

Q 3 = -{2(l+g )*o 

+ 2(j + q )s 

- (2j + Aq 2 + I2q + 6) Jo 

- (30gg + 84gg + 78®, + 24)} , (32) 

R = (jo - go - 2)Z 

- (3s - 2j + 6gg + 50g + 40)s 

+ [(3<fo + 10)j + llgg + Aqo - 18] j 

- (3gg + Mql + 246go + 104) . (33) 
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Eqs. (|23f - (133)) make it possible to estimate the present 
day values of f{R) and its first three derivatives as func- 
tion of the Hubble constant Hq and the cosmographic 
parameters (qo, jo, sq, Iq) provided a value for the mat- 
ter density parameter is given. This is a somewhat 
problematic point. Indeed, while the cosmographic pa- 
rameters may be estimated in a model independent way, 
the fiducial value for VLjm is usually the outcome of fit- 
ting a given dataset in the framework of an assumed dark 
energy scenario. However, it is worth noting that differ- 
ent models all converge towards the concordance value 
f2jvf ~ 0.25 which is also in agreement with astrophys- 
ical (model independent) estimates from the gas mass 
fraction in galaxy clusters. On the other hand, it has 
been proposed that f(R) theories may avoid the need for 
dark matter in galaxies and galaxy clusters [34| . In such 
a case, the total matter content of the universe is essen- 
tially equal to the baryonic one. According to the primor- 
dial elements abundance and the standard BBN scenario, 
we therefore get VIm — 0Jb/h 2 with uif, = fl^h 2 ~ 0.0214 
[35| and h the Hubble constant in units of lOOkm/s/Mpc. 
Setting h = 0.72 in agreement with the results of the HST 
Key project [3(|, we thus get Oj^ = 0.041 for a baryons 
only universe. We will therefore consider in the following 
both cases when numerical estimates are needed. 

It is worth noticing that Hq only plays the role of 
a scaling parameter giving the correct physical dimen- 
sions to f{R) and its derivatives. As such, it is not 
surprising that we need four cosmographic parameters, 
namely (qo, jo, s o,h), to fix the four f(R) related quan- 
tities f(R ), f'(Ro), f"(Ro), f"'(Ro)- It is also worth 
stressing that Eqs. (f2"3"|) - (|2^|) are linear in the f(R) quan- 
tities so that (qo,jo,so,lo) uniquely determine the for- 
mer ones. On the contrary, inverting them to get the 
cosmographic parameters as function of the f(R) ones, 
we do not get linear relations. Indeed, the field equa- 
tions in f(R) theories are nonlinear fourth order differ- 
ential equations in the scale factor a(t) so that fixing the 
derivatives of f(R) up to third order makes it possible 
to find out a class of solutions, not a single one. Each 
one of these solutions will be characterized by a different 
set of cosmographic parameters thus explaining why the 
inversion of Eqs. (|2"3"|) - (f3"3"|) does not give a unique result 
for (q ,jo,s ,lo)- 

As a final comment, we reconsider the underlying as- 
sumptions leading to the above derived relations. While 
Eqs. (fl7|) and (fl"8|) are exact relations deriving from a rig- 
orous application of the field equations, Eq. (f2"2")l heavily 
relies on having approximated f(R) with its third order 
Taylor expansion (|2"Tj) . If this assumption fails, the sys- 
tem should not be closed since a fifth unknown parameter 
enters the game, namely f( m '(Ro). Actually, replacing 
f(R) with its Taylor expansion is not possible for all class 
of f(R) theories. As such, the above results only hold in 
those cases where such an expansion is possible. More- 
over, by truncating the expansion to the third order, we 
are implicitly assuming that higher order terms are neg- 
ligible over the redshift range probed by the data. That 



is to say, we are assuming that : 



f^(Ro)(R-Ro) n « V J - { -P±(R- R ) m forn > 4 

^— ' ml 

(34) 

over the redshift range probed by the data. Checking 
the validity of this assumption is not possible without 
explicitly solving the field equations, but we can guess 
an order of magnitude estimate considering that, for all 
viable models, the background dynamics should not differ 
too much from the ACDM one at least up to z ~ 2. Us- 
ing then the expression of H(z) for the ACDM model, it 
is easily to see that R/Rq is a quickly increasing function 
of the redshift so that, in order Eq. (f34|) holds, we have to 
assume that / (,l) (^o) << f"'(Ro) for n > 4. This condi- 
tion is easier to check for many analytical f(R) models. 

Once such a relation is verified, we have still to worry 
about Eq. (|19p relying on the assumption that the cos- 
mological gravitational constant is exactly the same as 
the local one, i.e. the same as the one measured in the 
laboratory and entering the Newtonian Poisson equation. 
Actually, the cosmological gravitational constant should 
be identified with the one entering the perturbation equa- 
tions for a given f(R) model. Comparing the Newtonian 
Gjv and this cosmological G, one could infer whether the 
G entering the background equations is the same as the 
local one. Although this is outside our aims here, we can, 
in a first reasonable approximation, argue that the con- 
dition Gi OC ai = G cosmo could be replaced by the weaker 
relation G e ff(z = 0) = G(l + e) with e << 1. In this 
case, we should repeat the derivation of Eqs. (f2"3")) - (l2l?l) 
now using the condition f'(Ro) = (1 + e) _1 - Taylor ex- 
panding the results in e to the first order and comparing 
with the above derived equations, we can estimate the 
error induced by our assumption e = 0. The resulting 
expressions are too lengthy to be reported and depend in 
a complicated way on the values of the matter density pa- 
rameter Qm, the cosmographic parameters (qo,jo, sq,Iq) 
and e. However, we have numerically checked that the 
error induced on f(Ro), f"(Ro), f"'(Ro) are much lower 
than 10% for value of e as high as an unrealistic e ~ 0.1. 
We are therefore confident that our results are reliable 
also under such conditions. 



V. f(R) DERIVATIVES AND CPL MODELS 

A determination of f(R) and its derivatives in terms 
of the cosmographic parameters need for an estimate of 
these latter from the data in a model independent way. 
Unfortunately, even in the nowadays era of precision cos- 
mology, such a program is still too ambitious to give use- 
ful constraints on the f(R) derivatives, as we will see 
later. On the other hand, the cosmographic parameters 
may also be expressed in terms of the dark energy den- 
sity and EoS parameters so that we can work out what 
are the present day values of f(R) and its derivatives 
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giving the same (qo,jo,so,lo) of the given dark energy 
model. To this aim, it is convenient to adopt a param- 
eterized expression for the dark energy EoS in order to 
reduce the dependence of the results on any underlying 
theoretical scenario. Following the prescription of the 
Dark Energy Task Force [37| , we will use the Chevallier - 
Polarski- Linder (CPL) parameterization for the EoS set- 
ting [Hi : 

w = Wq + w a (l — a) = wo + w a z(l + z)^ 1 (35) 

so that, in a flat universe filled by dust matter and 
dark energy, the dimensionless Hubble parameter E(z) — 
H/Hq reads: 



the present day values of f(R) and its first three deriva- 
tives as function of (£Im> wo,w a ). It is worth noting that 
the f(R) model thus obtained is not dynamically equiv- 
alent to the starting CPL one. Indeed, the two models 
have the same cosmographic parameters only today. As 
such, for instance, the scale factor is the same between 
the two theories only over the time period during which 
the fifth order Taylor expansion is a good approxima- 
tion of the actual a(t). It is also worth stressing that 
such a procedure does not select a unique f(R) model, 
but rather a class of fourth order theories all sharing the 
same third order Taylor expansion of f{R). 

A. The AC DM case 



E 2 (z) = fi A f(l + 2) 3 + fix(l + ^) 3(1+ "' n+t, ' a) f 



(36) 



with Qx = 1 — f^M because of the flatness assumption. 
In order to determine the cosmographic parameters for 
such a model, we avoid integrating H(z) to get a(t) by 
noting that d/dt = —(1 + z)H(z)d/dz. We can use such 

a relation to evaluate (H, H, H, d A H/dt A ) and then solve 
Eqs.([2])- ©, evaluated in z = 0, with respect to the pa- 
rameters of interest. Some algebra finally gives : 



1 3 



(37) 



With these caveats in mind, it is worth considering first 
the ACDM model which is obtained by setting (wq, w a ) — 
(—1,0) in the above expressions thus giving: 



<7o = 

jo = 

so = 

la = 



1- 2 fi M 



27 

1 + 30 M + — n 2 M 



(41) 



Jo 



-(1 - n M ) [3W (1 + Wq) + W a 



(38) 



When inserted into the expressions for the f(R) quanti- 
ties, these relations give the remarkable result : 



so 



7 33 n o ^ 

~2 - -ji 1 - ^M)W a 
9 

-(i - n M ) [9 + (7 - n M ) Wa ] w 



4 
27 



(l-fiA,)(16-3fi M K 



(1 - Q M )(3 - n M )w 3 , 



h = y + LJlK [213 + (7 - Sl M )w a ] w a 



1 - 



M ) 



9(82 - 2m M )w a ) wo 



f(Ro) =R a + 2A , f"(R a ) = f"'(Ro) = , (42) 

so that we obviously conclude that the only f(R) theory 
having exactly the same cosmographic parameters as the 
ACDM model is just f(R) oc R, i.e. GR. It is worth 
noticing that such a result comes out as a consequence 
of the values of (qo,jo) m the ACDM model. Indeed, 
(39) should we have left (sq,Iq) undetermined and only fixed 
(<Zo> jo) to the values in (|4"Tj) . we should have got the same 
result in fl2|. Since the ACDM model fits well a large 
set of different data, we do expect that the actual values 
of {qo,jo, Sq, lo) do not differ too much from the ACDM 
ones. Therefore, we plug into Eqs. ([2"3"]) - (|33|) the following 
expressions : 



67-21fijtf + -(23-llQ M H 



27 

+ — (l-n M )(47-24fi M )^ 
81 

+ — (i - o M )(3 - m M )wt . 



(40) 



Inserting Eqa.((57 |> - P? | into Eqs.([2"3" ]) - ([55 ]) . we get 
lengthy expressions (which we do not report here) giving 



q a = qox(l + e q ) , j =jo x ( 1 + £ i) 



s = Sox{l + e s ) , l Q = lo x(l +e/) , 

with (qo,jo, s oJo) S iven b y Eqs.flU]) and (e q , Ej, e s , £;) 
quantifyin the deviations from the ACDM values allowed 
by the data. A numerical estimate of these quantities 
may be obtained, e.g., from a Markov chain analysis, 
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but this is outside our aims. Since we are here inter- 
ested in a theoretical examination, we prefer to consider 
an idealized situation where the four quantities above 
all share the same value e << 1. In such a case, we 
can easily investigate how much the corresponding f(R) 
deviates from the GR one considering the two ratios 
f"(Ro)/f(Ro) and f'"(R )/f(R ). Inserting the above 
expressions for the cosmographic parameters into the 
exact (not reported) formulae for f(Ro), f"(Ro) and 
f"'(Ro), taking their ratios and then expanding to first 
order in e, we finally get : 



V20 



mo 



64 - 60 M (90 M + ; 



[3(90 M + 74)0 M - 556] 0| 



M 



X — 

16 27 



(43) 



6 [(810m - 110)O M + 40] n M + 16 
[3(9fi M + 74)fi 



M 



556] n 2 M + 16 



2430 M 



(44) 

having defined ij 20 = f"(Ro)/f(R )xH$ and 7730 = 
f"'(Ro)/f(Ro)xHQ which, being dimensionless quanti- 
ties, are more suited to estimate the order of magnitudes 
of the different terms. Inserting our fiducial values for 
fl M , we get: 

7720 ^ 0.15 x s for SIm = 0.041 
7/20 ^ -0.12 x e for fl M = 0.250 



mo 



mo 



4 x 



for fi 



M 



0.041 



-0.18 x e for 9, M = 0.250 



For values of e up to 0.1, the above relations show that 
the second and third derivatives are at most two orders 
of magnitude smaller than the zeroth order term f(Rg). 
Actually, the values of 7730 for a baryon only model (first 
row) seems to argue in favor of a larger importance of the 
third order term. However, we have numerically checked 
that the above relations approximates very well the exact 
expressions up to e ~ 0.1 with an accuracy depending on 
the value of CIm, being smaller for smaller matter density 
parameters. Using the exact expressions for 7720 and 7730, 
our conclusion on the negligible effect of the second and 
third order derivatives are significantly strengthened. 

Such a result holds under the hypotheses that the nar- 
rower are the constraints on the validity of the ACDM 
model, the smaller are the deviations of the cosmographic 
parameters from the ACDM ones. It is possible to show 
that this indeed the case for the CPL parametrization 
we are considering. On the other hand, we have also as- 
sumed that the deviations (eg, Sj, s s , £;) take the same 
values. Although such hypothesis is somewhat ad hoc, 
we argue that the main results are not affected by giving 
it away. Indeed, although different from each other, we 
can still assume that all of them are very small so that 
Taylor expanding to the first order should lead to addi- 
tional terms into Eqs. (|4"3"l) - (I44[) which are likely of the 



same order of magnitude. We may therefore conclude 
that, if the observations confirm that the values of the 
cosmographic parameters agree within ~ 10% with those 
predicted for the ACDM model, we must conclude that 
the deviations of f(R) from the GR case, f(R) oc R, 
should be vanishingly small. 

It is worth stressing, however, that such a conclusion 
only holds for those f(R) models satisfying the constraint 
(|34)) . It is indeed possible to work out a model having 
f(Ro) oc Ro, f'(R Q ) = f"(R ) = 0, but f^(R ) # for 
some n. For such a (somewhat ad hoc) model, Eq. (j3~4"|) 
is clearly not satisfied so that the cosmographic param- 
eters have to be evaluated from the solution of the field 
equations. For such a model, the conclusion above does 
not hold so that one cannot exclude that the resulting 
(<Zo, jo, so, lo) are within 10% of the ACDM ones. 



B. The constant EoS model 

Let us now take into account the condition w = — 1, 
but still retains w a = thus obtaining the so called 
quiessence models. In such a case, some problems arise 
because both the terms (jo — <7o — 2) and 1Z may van- 
ish for some combinations of the two model parameters 
(Omj wo). For instance, we find that jo — qo — 2 = for 
Wo = (wi, W2) with : 



1 



Wl 



w 2 



1 - o M + A /(i-fiM)(4-n M ) 



4-0 



M 



V(1-O m )(4-O m ) 



On the other hand, the equation TZ(Qm , wq) = may 
have different real roots for w depending on the adopted 
value of Am. Denoting collectively with w nu ii the values 
of wo that, for a given Om, make (jo — qo — 2)7?.(Om, Wo) 
taking the null value, we individuate a set of quiessence 
models whose cosmographic parameters give rise to di- 
vergent values of f(Ro, f"(Ro) and f"'(Ro), For such 
models, f(R) is clearly not defined so that we have to 
exclude these cases from further consideration. We only 
note that it is still possible to work out a f(R) theory re- 
producing the same background dynamics of such mod- 
els, but a different route has to be used. 

Since both qo and jo now deviate from the ACDM val- 
ues, it is not surprising that both f"(Ro) and f'"(Ro) 
take finite non null values. However, it is more interest- 
ing to study the two quantities 7720 and 7730 defined above 
to investigate the deviations of f(R) from the GR case. 
These are plotted in Figs.[T] and [5] for the two fiducial 
Om values. Note that the range of wo in these plots have 
been chosen in order to avoid divergences, but the lessons 
we will draw also hold for the other wo values. 

As a general comment, it is clear that, even in this case, 
f"(Ro) and f"'(Ro) are from two to three orders of mag- 
nitude smaller that the zeroth order term f(Ro). Such a 



-0.015 - x 
-0.02 



-1.3 -1.2 -1.1 -1 

wo 



-0.9 -0.8 -0.7 



FIG. 1: The dimensionless ratio between the present day val- 
ues of f"(R) and f(R) as function of the constant EoS wo of 
the corresponding quiessence model. Short dashed and solid 
lines refer to models with Qm = 0.041 and 0.250 respectively. 



FIG. 2: The dimensionless ratio between the present day val- 
ues of f"'(R) and f(R) as function of the constant EoS too of 
the corresponding quiessence model. Short dashed and solid 
lines refer to models with Qm = 0.041 and 0.250 respectively. 



result could be yet guessed from the previous discussion 
for the ACDM case. Actually, relaxing the hypothesis 
wq = —1 is the same as allowing the cosmographic pa- 
rameters to deviate from the ACDM values. Although 
a direct mapping between the two cases cannot be es- 
tablished, it is nonetheless evident that such a relation 
can be argued thus making the outcome of the above 
plots not fully surprising. It is nevertheless worth not- 
ing that, while in the ACDM case, 7720 and 7730 always 
have opposite signs, this is not the case for quiessence 
models with w > — 1. Indeed, depending on the value of 
CImi we can have f(R) theories with both 7720 and 7730 
positive. Moreover, the lower is U~m, the higher are the 
ratios 7720 and 7730 for a given value of wq. This can be 
explained qualitatively noticing that, for a lower SIm, the 
density parameter of the curvature fluid (playing the role 
of an effective dark energy) must be larger thus claiming 
for higher values of the second and third derivatives (see 
also [39| for a different approach to the problem). 



C. The general case 

Finally, we consider evolving dark energy models with 
w a ^ 0. Needless to say, varying three parameters allows 
to get a wide range of models that cannot be discussed 
in detail. Therefore, we only concentrate on evolving 
dark energy models with wo = —1 in agreement with 
some most recent analysis. The results on 7720 and 7730 
are plotted in Figs. [3] and Q] where these quantities as 
functions of w a . Note that we are considering models 
with positive w a so that w(z) tends to wq + w a > wq 
for z — > 00 so that the EoS dark energy can eventually 
approach the dust value w = 0. Actually, this is also the 
range favored by the data. We have, however, excluded 
values where 7720 or 7730 diverge. Considering how they 
are defined, it is clear that these two quantities diverge 
when f(Ro) — so that the values of (wo,w a ) making 



(f?20 , 7730 ) to diverge may be found solving : 
V {w 0l w a )Q M + Qo(w 0l w a ) = 

where Vo(wo, w a ) and Qo(wo, w a ) are obtained by insert- 
ing Eqs.(|37l) - (J40J) into the defintions (27]) - ([28]). For such 
CPL models, there is no any f(R) model having the same 
cosmographic parameters and, at the same time, satisfy- 
ing all the criteria needed for the validity of our proce- 
dure. Actually, if /(i?o) = 0, the condition is likely 
to be violated so that higher than third order must be in- 
cluded in the Taylor expansion of f(R) thus invalidating 
the derivation of Eqs.(l3]) - ([26]) . 

Under these caveats, Figs.|3]and[2]demonstrate that al- 
lowing the dark energy EoS to evolve does not change sig- 
nificantly our conclusions. Indeed, the second and third 
derivatives, although being not null, are nevertheless neg- 
ligible with respect to the zeroth order term thus argu- 
ing in favour of a GR-like f(R) with only very small 
corrections. Such a result is, however, not fully unex- 
pected. From Eqs. ([37|) and (f38|) . we see that, having 
setted wq — —1, the qq parameter is the same as for the 
ACDM model, while j reads jfr + (3/2) (1 - Sl M )ic„. As 
we have stressed above, the Hilbert - Einstein Lagrangian 
f(R) = R+ 2 A is recovered when (go, jo) = (lodo ) what- 
ever the values of (sq, lo) are. Introducing a!o o /0 makes 
(so, lo) to differ from the ACDM values, but the first two 
cosmographic parameters are only mildly affected. Such 
deviations are then partially washed out by the compli- 
cated way they enter in the determination of the present 
day values of f(R) and its first three derivatives. 

VI. CONSTRAINING f(R) PARAMETERS 

In the previous section, we have worked an alternative 
method to estimate f(Ro), f"(Ro), f"'(Ro) resorting to 
a model independent parameterization of the dark en- 
ergy EoS. However, in the ideal case, the cosmographic 
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FIG. 3: The dimensionless ratio between the present day val- 
ues of f" (Ft) and f(R) as function of the w a parameter for 
models with wo — — 1. Short dashed and solid lines refer to 
models with SIjvj = 0.041 and 0.250 respectively. 



FIG. 4: The dimensionless ratio between the present day val- 
ues of f'" (Ft) and f(Ft) as function of the w a parameter for 
models with wo = — 1. Short dashed and solid lines refer to 
models with Qm — 0.041 and 0.250 respectively. 



parameters are directly estimated from the data so that 
Eqs. (j2"3")) - ([55]) can be used to infer the values of the f(R) 
related quantities. These latter can then be used to put 
constraints on the parameters entering an assumed fourth 
order theory assigned by a f(R) function characterized by 
a set of parameters p = (pi, . . . ,p n ) provided that the hy- 
potheses underlying the derivation of Eqs. ([2"3")) - are 
indeed satisfied. We show below two interesting cases 
which clearly highlight the potentiality and the limita- 
tions of such an analysis. 



A. Double power law Lagrangian 

As a first interesting example, we set : 



f(R) =R(l + aR n + 0R~ 



(45) 



with n and m two positive real numbers (see, for exam- 
ple, [13] for some physical motivations). The following 
expressions are immediately obtained : 



( f(R ) = Ro(l + aR% + 0R^> 



f'(Ro) = 1 + a (n + 1)R$ - 0(m - 1)R, 



— m 




(1+m) 



f"(Ro) = an(n + 1)-R£ _1 + 0m(m - 1)R~ 

f'"(Ro) = an(n+l)(n-l)R%- 2 

- f3m(m + l)(m- l)R^ (2+m) 



Denoting by <pi (with i = 0, . . . ,3) the values of f^(Ro) 
determined through Eqs. (f23|) - ([33]) . we can solve : 



f(Ro) = ^ 

f'(Ra) = 0i 

f"(R ) = 02 

f"'(Ro) = 03 



which is a system of four equations in the four unknowns 
(a, (3, n, m) that can be analytically solved proceeding as 
follows. First, we solve the first and second equation with 
respect to (a, [3) obtaining : 



n + m \ Ko 



(46) 







1 + n 



1 



7?" 



n + m \ Rq 
while, solving the third and fourth equations, we get : 



= 



<hK~ n [l + m+(<fe/0 2 )flo] 
n(n + l)(n + m) 

<j> 2 Rl +n [l-n+ {</> 3 /<h)Ro] 
m(l — m)(n + m) 



(47) 



Equating the two solutions, we get a systems of two equa- 
tions in the two unknowns (n, m), namely: 



n(n + l)(l-m) (1-fa/Rg) 
<foRa [l + m+ (h/h)Ro\ 

m(n + l)(m-l)(l-0 o /i?o) 
<f> 2 Ro [l-n+(fa/<h)Ro] 



= 1 



= 1 



(48) 



Solving with respect to m, we get two solutions, the first 
one being m — —n which has to be discarded since makes 
(a, 0) goes to infinity. The only acceptable solution is : 



m 



[l-n+ (fo/fa)Ro] 



(49) 



which, inserted back into the above system, leads to a 
second order polynomial equation for n with solutions : 
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with i?oc = Ro I ' Rc and : 



1 

where we have defined : 



b 2 O± 2 i?o(l 



00 /Rq 



(50) 



+ 6 (i? o 0o - 2i?o^o + Rl) 
+ 9 (0g - 2i? o 0o + i?o) 02 

+ 4 (i?^ - fig) 03 . 



"2^3 



(51) 



Depending on the values of (<?o, jo, so, ^o)> Eg. ([50")) may 
lead to one, two or any acceptable solution, i.e. real 
positive values of n. This solution has then to be in- 
serted back into Eq. (|4"9"|) to determine m and then into 
Eqs. (|4U)) or (ITT)) to estimate (a,/?). If the final values of 
(a, /3, n, m) are physically viable, we can conclude that 
the model in Eq. (|4"5")) is in agreement with the data giv- 
ing the same cosmographic parameters inferred from the 
data themselves. Exploring analytically what is the re- 
gion of the ((/o, jo, sq,Io) parameter space which leads to 
acceptable (a, (3, n, m) solutions is a daunting task far 
outside the aim of the present paper. 



B. HS model 

One of the most pressing problems of f(R) theories is 
the need to escape the severe constraints imposed by the 
Solar System tests. A successful model has been recently 
proposed by Hu and Sawicki [2l[ (HS) setting 8 : 



-f3 2 R& + 4(3R%c 
3(1 - 0*F&) 

-2(1 - mz c ) 2 



(54) 



Equating Eqs. (f53")) to the four quantities (0o, 0i, 02, 03) 
defined as above, we could, in principle, solve this system 
of four equations in four unknowns to get (a, 0, R c , n) in 
terms of (0o, 01, 02, 03) and then, using Eqs. (f2"3"l) - (|3"3"j) as 
functions of the cosmographic parameters. However, set- 
ting 0i = 1 as required by Eq. (|24p gives the only trivial 
solution anR c = so that the HS model reduces to the 
Einstein- Hilbcrt Lagrangian f(R) = R. In order to es- 
cape this problem, we can relax the condition /'(i?o) = 1 
to f'(Ro) = (1 + e) -1 - As we have discussed in Sect. IV, 
this is the same as assuming that the present day ef- 
fective gravitational constant G e ff,o = Gn/ f'(Ro) only 
slightly differs from the usual Newtonian one which seems 
to be a quite reasonable assumption. Under this hypoth- 
esis, we can analytically solve for (a, (3, R c , n) in terms of 
(0o, e, 02, 03). The actual values of (0o, 02, 03) will be no 
more given by Eqs. (|2"3")l - (j2l)|) , but we have checked that 



they deviate from those expressions 9 much less than 10% 
for e up to 10% well below any realistic expectation. 
With this caveat in mind, we first solve 



f(Ro) = <h , /"(iJo) = (l + e)- 

to get : 



f(R) = R-R C 



a{R/R c ) n 
1 + (3(R/R c ) n 



(52) 



As for the double power law model discussed above, there 
are four parameters which we can be expressed in terms 
of the cosmographic parameters (go, jo, so, lo)- 
As a first step, it is trivial to get : 



f(Ro) 
f(Ro) 
f'(Ro) 
f"(Ro) 



r r aR °c 



1 + Wc 

_ anR c R^ c 
Rq(1 + {3R% c ) 2 

anRcRfic [(1 - n) + j3(l + n)R a l c 
R 2 o(l + (3RS c ) 3 

anR c R^ c (An 2 +Bn + C) 

RW+mJ 4 



(53) 



/3 = 



n{l+e) 

e 

n(l +e) 



Rq 
Rc 

Rq 
Rc 



Ro 



1 _ x?_ 



Ro n(l + e) 



Inserting these expressions in Eqs. (|55|) . it is easy to check 
that R c cancels out so that we can no more determine its 
value. Such a result is, however, not unexpected. Indeed, 
Eq. ([52]) can trivially be rewritten as : 



f(R) = R 



&R r ' 



l + f3R n 



with q = aRl~ n and (3 = (3R~ n which are indeed the 
quantities that are determined by the above expressions 
for (a, /3). Reversing the discussion, the present day val- 
ues of /W (_R) depend on (a, (3, R c ) only through the two 



Note that such a model does not pass the matter instability test 
so that some viable generalizations [SI have been proposed. 



Note that the correct expressions for (phio , ip2 , 4*3 ) ma y still for- 
mally be written as Eqs. (12311 - (1261) . but the polynomials entering 
them are now different and also depend on powers of e. 
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parameters (a, (3). As such, the use of cosmographic pa- 
rameters is unable to break this degeneracy. However, 
since R c only plays the role of a scaling parameter, we 
can arbitrarily set its value without loss of generality. 

On the other hand, this degeneracy allows us to get 
a consistency relation to immediately check whether the 
HS model is viable or not. Indeed, solving the equation 
f"(Ro) = fa, we get: 

ifa/Ro) + [(1 + e)/e](l - faRo) - (1 - e)/(l + e) 

n = ; 7— , 

1 - fa/R 

which can then be inserted into the equations f"'(Ro) = 
fa to obtain a complicated relation among {fa, fa, ^3) 
which we do not report for sake of shortness. Solving such 
a relation with respect to fa / fa and Taylor expanding to 
first order in e, the constraint we get reads : 



1 



Ro 



1 



1 



2s 



1 - fa/Ro 



If the cosmographic parameters {qo, jo, so,lo) are known 
with sufficient accuracy, one could compute the values of 
{Ro, fa, fa-fa) for a given e (eventually using the expres- 
sions obtained for e = 0) and then check if they satisfied 
this relation. If this is not the case, one can immediately 
give off the HS model also without the need of solving the 
field equations and fitting the data. Actually, given the 
still large errors on the cosmographic parameters, such a 
test only remains in the realm of (quite distant) future 
applications. However, the HS model works for other 
tests as shown in [2l| and so a consistent cosmography 
analysis has to be combined with them. 



where CV, are the elements of the covariance matrix (be- 
ing C u = a*.), we have set {pi,P2,P3,Pa) = {Qo,3o,so,Io)- 
and assumed that the erorr <jm on Qm is uncorrelated 
with those on p. Note that this latter assumption strictly 
holds if the matter density parameter is estimated from 
an astrophysical method (such as estimating the total 
matter in the universe from the estimated halo mass 
function). Alternatively, we will assume that Om is 
constrained by the CMBR related experiments. Since 
these latter mainly probes the very high redshift uni- 
verse {z ~ ziss ~ 1089), while the cosmographic pa- 
rameters are concerned with the present day cosmo, one 
can argue that the determination of Qm is not affected 
by the details of the model adopted for describing the 
late universe. Indeed, we can reasonably assume that, 
whatever is the dark energy candidate or f(R) theory, 
the CMBR era is well approximated by the standard GR 
with a model comprising only dust matter. As such, we 
will make the simplifying (but well motivated) assump- 
tion that o m may be reduced to very small values and is 
uncorrelated with the cosmographic parameters. 

Under this assumption, the problem of estimating the 
errors on g(f2jy,p) reduces to estimating the covariance 
matrix for the cosmographic parameters given the details 
of the data set used as observational constraints. We 
address this issue by computing the Fisher information 
matrix (see, e.g., [42j and references therein) defined as: 



F, 



d 2 L 



(56) 



VII. CONSTRAINTS ON f(R) DERIVATIVES 
FROM THE DATA 

Eqs. (f2"3"l) - ([55)1 relate the present day values of f{R) 
and its first three derivatives to the cosmographic pa- 
rameters (<?o, jo, so,lo) and the matter density £Im- In 
principle, therefore, a measurement of these latter quan- 
tities makes it possible to put constraints on /w(J2o), 
with i = {0, . . . , 3}, and hence on the parameters of a 
given fourth order theory through the method shown in 
the previous section. Actually, the cosmographic param- 
eters are affected by errors which obviously propagate 
onto the f{R) quantities. Actually, the covariance ma- 
trix for the cosmographic parameters is not diagonal so 
that one has also take care of this to estimate the final 
errors on 

/W(i?o)- A similar discussion also holds for the 
errors on the dimcnsionless ratios 7720 and 7730 introduced 
above. As a general rule, indicating with g{QM,p) a 
generic f(R) related quantity depending on Qm and the 
set of cosmographic parameters p, its uncertainty reads : 



dg 



on 



A I 



M 



i=4 

E 



dg 



dpi 



Y~2 — —C- 

ffi dpi dpj 13 



withL = -2 In £(0i, . . . ,6 n ), £(6>i,.. .,0 n ) the likelihood 
of the experiment, (#1, . . . , 9 n ) the set of parameters to 
be constrained, and (. . .} denotes the expectation value. 
Actually, the expectation value is computed by evaluat- 
ing the Fisher matrix elements for fiducial values of the 
model parameters (#1, . . . , 6 n ), while the covariance ma- 
trix C is finally obtained as the inverse of F. 

A key ingredient in the computation of F is the defini- 
tion of the likelihood which depends, of course, of what 
experimental constraint one is using. To this aim, it is 
worth remembering that our analysis is based on fifth 
order Taylor expansion of the scale factor a(t) so that 
we can only rely on observational tests probing quan- 
tities that are well described by this truncated series. 
Moreover, since we do not assume any particular model, 
we can only characterize the background evolution of the 
universe, but not its dynamics which, being related to the 
evolution of perturbations, unavoidably need the specifi- 
cation of a physical model. As a result, the SNela Hubble 
diagram is the ideal test 10 to constrain the cosmographic 
parameters. We therefore defined the likelihood as : 



(55) 



See the conclusions for further discussion on this issue. 
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£(H ,p) cxexp- x 2 (^o,p)/2 



X 2 (ffo,p)=E^r° 



where the distance modulus to redshift z reads : 



(57) 



Hth(z,H ,p) = 25 + 51og(c/H ) + 51ogd L (z,p) , (58) 
and cIl(z) is the Hubble free luminosity distance : 



d L (z) = (l + z) 



dz 



o H(z)/H 



(59) 



Using the fifth order Taylor expansion of the scale factor, 
we get for di,(z,p) an analytical expression (reported in 
Appendix A) so that the computation of does not 
need any numerical integration (which makes the esti- 
mate faster). As a last ingredient, we need to specify the 
details of the SNela survey giving the redshift distribu- 
tion of the sample and the error on each measurement. 
Following [43(| , we adopt 11 : 



a ( z) 



with Zmax the maximum redshift of the survey, cr sys an 
irreducible scatter in the SNela distance modulus and a m 
to be assigned depending on the photometric accuracy. 

In order to run the Fisher matrix calculation, we have 
to set a fiducial model which we set according to the 
ACDM predictions for the cosmographic parameters. For 
£Im = 0.3 and h = 0.72 (with h the Hubble constant in 
units of lOOkm/s/Mpc), we get: 

(«b,Jo, *o,*o) = (-0.55,1.0,-0.35,3.11) . 

As a first consistency check, we compute the Fisher ma- 
trix for a survey mimicking the recent database in Q 
thus setting (AfsNeia,cr m ) = (192,0.33). After marginal- 
izing over h (which, as well known, is fully degenerate 
with the SNela absolute magnitude M), we get for the 
uncertainties : 

(ai, ©a, cr 3 , <j 4 ) = (0.38, 5.4, 28.1, 74.0) 

where we are still using the indexing introduced above 
for the cosmographic parameters. These values compare 



11 Note that, in [431 . the authors assume the data are separated 
in redshift bins so that the error becomes 

Nbi n (z I Zmax) 2 &m with Mbin the number of SNela in a bin. How- 
ever, we prefer to not bin the data so that Mbin = !• 



reasonably well with those obtained from a cosmographic 
fitting of the Gold SNela dataset 12 : 



3o 



-0.90±0.65 , j = 2.7±6.7 



s = 36.5±52.9 , l = 142.7±320 . 

Because of the Cramer - Rao theorem, the Fisher matrix 
approach is known to provide the minimum variance er- 
rors a given experiment can attain thus giving higher 
limits to its accuracy on the determination of a set of 
parameters. This is indeed the case with the comparison 
suggesting that our predictions are quite optimistic. It is 
worth stressing, however, that the analysis in [44| used 
the Gold dataset which is poorer in high z SNela than 
the 0] one we are mimicking so that larger errors on the 
higher order parameters (soj^o) are expected. 

Rather than computing the errors on f(Ro) and its 
first three derivatives, it is more interesting to look at the 
precision attainable on the dimcnsionless ratios (??20)^30 
introduced above since they quantify how much devia- 
tions from the linear order are present. For the fiducial 
model we are considering, both 7720 and 7730 vanish, while, 
using the covariance matrix for a present day survey and 
setting ctm/^m — 10%, their uncertainties read: 

(<720,CT 3 o) - (0.04,0.04) . 

As an application, we can look at Figs.Q] and [2] show- 
ing how (?720j»73o) depend on the present day EoS wq for 
f(R) models sharing the same cosmographic parameters 
of a dark energy model with constant EoS. As it is clear, 
also considering only the la range, the full region plot- 
ted is allowed by such large constraints on (7720, %o) thus 
meaning that the full class of corresponding f(R) theo- 
ries is viable. As a consequence, we may conclude that 
the present day SNela data are unable to discriminate 
between a A dominated universe and this class of fourth 
order gravity theories. 

As a next step, we consider a SNAP -like survey [45] ] 
thus setting {M S Neia,o m ) = (2000,0.02). We use the 
same redshift distribution in Table 1 of [43[ and add 
300 nearby SNela in the redshift range (0.03,0.08). The 
Fisher matrix calculation gives for the uncertainties on 
the cosmographic parameters : 

((71,0-2, 0-3, 0-4) = (0.08,1.0,4.8,13.7) . 



12 Actually, such estimates have been obtained computing the mean 
and the standard deviation from the marginalized likelihoods of 
the cosmographic parameters. As such, the central values do not 
represent exactly the best fit model, while the standard devia- 
tions do not give a rigorous description of the error because the 
marginalized likelihoods are manifestly non - Gaussian. Never- 
theless, we are mainly interested in an order of magnitude esti- 
mate so that we do not care about such statistical details. 
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The significant improvement of the accuracy in the de- 
termination of (go> JO) s oJo) translates in a reduction of 
the errors on (7720,7730) which now read: 

(<T20, ^3o) = (0.007,0.008) 

having assumed that, when SNAP data will be available, 
the matter density parameter Qm has been determined 
with a precision ctm/^m ~ 1%- Looking again at Figs. [I] 
and [21 it is clear that the situation is improved. Indeed, 
the constraints on 7720 makes it possible to narrow the 
range of allowed models with low matter content (the 
dashed line) , while models with typical values of fl m are 
still viable for wo covering almost the full horizontal axis. 
On the other hand, the constraint on 7730 is still too weak 
so that almost the full region plotted is allowed. 

Finally, we consider an hypothetical future SNela 
survey working at the same photometric accuracy as 
SNAP and with the same redshift distribution, but in- 
creasing the number of SNela up to MsNeia = 6xl0 4 
as expected froim_ e.g., DES 0, PanSTARRS [47[, 
SKYMAPPER [Hi, while still larger numbers may po- 
tentially be achieved by ALPACA H| and LSST [5cj. 
Such a survey can achieve : 

(<ri, OS, c-3, o*) = (0.02,0.2,0.9,2.7) 
so that, with ctm/Qm ~ 0.1%, we get : 

(^0,^30) = (0.0015,0.0016) . 

Fig.Q]shows that, with such a precision on 7720, the region 
of ttjo values allowed essentially reduces to the ACDM 
value, while, from Fig. [21 it is clear that the constraint on 
7730 definitively excludes models with low matter content 
further reducing the range of u>o values to quite small 
deviations from the wq = — 1 . We can therefore conclude 
that such a survey will be able to discriminate between 
the concordance ACDM model and all the f(R) theories 
giving the same cosmographic parameters as quiessence 
models other than the ACDM itself. 

A similar discussion may be repeated for f(R) models 
sharing the same (qo, jo, sq, Iq) values as the CPL model 
even if it is less intuitive to grasp the efficacy of the sur- 
vey being the parameter space multivalued. For the same 
reason, we have not explored what is the accuracy on the 
double power -law or HS models, even if this is techni- 
cally possible. Actually, one should first estimate the er- 
rors on the present day value of f(R) and its three time 
derivatives and then propagate them on the model pa- 
rameters using the expressions obtained in Sect. VI. The 
multiparameter space to be explored makes this exercise 
quite cumbersome so that we leave it for a furthcoming 
work where we will explore in detail how these models 
compare to the present and future data. 

VIII. CONCLUSIONS 

The recent amount of good quality data have given a 
new input to the observational cosmology. As often in 



science, new and better data lead to unexpected discov- 
eries as in the case of the nowadays accepted evidence for 
cosmic acceleration. However, a fierce and strong debate 
is still open on what this cosmic speed up implies for 
theoretical cosmology. The equally impressive amount 
of different (more or less) viable candidates have also 
generated a great confusion so that model independent 
analyses are welcome. A possible solution could come 
from the cosmography of the universe rather than as- 
suming ad hoc solutions of the cosmological Friedmann 
equations. Present day and future SNela surveys have 
renewed the interest in the determination of the cosmo- 
graphic parameters so that it is worth investigating how 
these quantities can constrain cosmological models. 

Motivated by this consideration, in the framework of 
metric formulation of f(R) gravity, we have here derived 
the expressions of the present day values of f(R) and its 
first three derivatives as function of the matter density 
parameter Om, the Hubble constant Hq and the cosmo- 
graphic parameters (go, jo, So, Iq). Although based on a 
third order Taylor expansion of f{R), we have shown 
that such relations hold for a quite large class of models 
so that they are valid tools to look for viable f(R) models 
without the need of solving the mathematically difficult 
nonlinear fourth order differential field equations. 

Notwithstanding the common claim that we live 
in the era of precision cosmology, the constraints on 
(<7o, Jo, so, Iq) are still too weak to efficiently apply the 
program we have outlined above. As such, we have shown 
how it is possible to establish a link between the popu- 
lar CPL parameterization of the dark energy equation 
of state and the derivatives of f(R), imposing that they 
share the same values of the cosmographic parameters. 
This analysis has lead to the quite interesting conclu- 
sion that the only f(R) function able to give the same 
values of (qo, jo, S0J0) as the ACDM model is indeed 
f(R) = R+2A. If future observations will tell us that the 
cosmographic parameters are those of the ACDM model, 
we can therefore rule out all f(R) theories satisfying the 
hypotheses underlying our derivation of Eqs. (|2"5)) - (|2"6"]) . 
Actually, such a result should not be considered as a no 
way out for higher order gravity. Indeed, one could still 
work out a model with null values of f"(Ro) and f'"(Ro) 
as required by the above constraints, but non - vanishing 
higher order derivatives. One could well argue that such 
a contrived model could be rejected on the basis of the 
Occam's razor, but nothing prevents from still taking it 
into account if it turns out to be both in agreement with 
the data and theoretically well founded. 

If new SNela surveys will determine the cosmographic 
parameters with good accuracy, acceptable constraints 
on the two dimensionless ratios 7720 oc f"(Ro)/f(Ro) and 
V30 oc f'"(Ro)/f{Ro) could be obtained thus allowing to 
discriminate among rival f(R) theories. To investigate 
whether such a program is feasible, we have pursued 
a Fisher matrix based forecasts of the accuracy future 
SNela surveys can achieve on the cosmographic parame- 
ters and hence on (7720,7730). It turns out that a SNAP- 
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like survey can start giving interesting (yet still weak) 
constraints allowing to reject f(R) models with low mat- 
ter content, while a definitive improvement is achievable 
with future SNela survey observing ~ 10 4 objects thus 
making it possible to discriminate between ACDM and a 
large class of fourth order theories. It is worth stressing, 
however, that the measurement of £Im should come out 
as the result of a model independent probe such as the 
gas mass fraction in galaxy clusters which, at present, is 
still far from the 1% requested precision. On the other 
hand, one can also rely on the Q,m estimate from the 
CMBR anisotropy and polarization spectra even if this 
comes to the price of assuming that the physics at recom- 
bination is strictly described by GR so that one has to 
limit its attention to f(R) models reducing to f(R) oc R 
during that epoch. However, such an assumption is quite 
common in many f(R) models available in literature so 
that it is not a too restrictive limitation. 

A further remark is in order concerning what kind 
of data can be used to constrain the cosmographic pa- 
rameters. The use of the fifth order Taylor expansion 
of the scale factor makes it possible to not specify any 
underlying physical model thus relying on the minimal- 
ist assumption that the universe is described by the flat 
Robertson - Walker metric. While useful from a theoreti- 
cal perspective, such a generality puts severe limitations 
to the dataset one can use. Actually, we can only resort 
to observational tests depending only on the background 
evolution so that the range of astrophysical probes re- 
duces to standard candles (such as SNela and possibly 
Gamma Ray Bursts) and standard rods (such as the an- 
gular size - redshift relation for compact radiosources). 
Moreover, pushing the Hubble diagram to z ~ 2 may 
rise the question of the impact of gravitational lensing 
amplification on the apparent magnitude of the adopted 
standard candle. The magnification probability distribu- 
tion function depends on the growth of perturbations [5l[ 
so that one should worry about the underlying physical 
model in order to estimate whether this effect biases the 
estimate of the cosmographic parameters. However, it 
has been shown [J, [52j that the gravitational lensing am- 
plification does not alter significantly the measured dis- 
tance modulus for z <~ 1 SNela. Although such an anal- 
ysis has been done for GR based models, we can argue 
that, whatever is the f(R) model, the growth of perturba- 
tions finally leads to a distribution of structures along the 
line of sight that is as similar as possible to the observed 
one so that the lensing amplification is approximately the 
same. We can therefore argue that the systematic error 
made by neglecting lensing magnification is lower than 
the statistical ones expected by the future SNela surveys. 
On the other hand, one can also try further reducing this 
possible bias using the method of flux averaging [53| even 
if, in such a case, our Fisher matrix calculation should 
be repeated accordingly. It is also worth noting that 
the constraints on the cosmographic parameters may be 
tigthened by imposing some physically motivated priors 
in the parameter space. For instance, we can impose that 



the Hubble parameter H(z) stays always positive over the 
full range probed by the data or that the transition from 
past deceleation to present acceleration takes place over 
the range probed by the data (so that we can detect it). 
Such priors should be included in the likelihood definition 
so that the Fisher matrix should be recomputed which is 
left for a forthcoming paper. 

Although the present day data are still too limited to 
efficiently discriminate among rival f(R) models, we are 
confident that an aggressive strategy aiming at a very 
precise determination of the cosmographic parameters 
could offer stringent constraints on higher order gravity 
without the need of solving the field equations or address- 
ing the complicated problems related to the growth of 
perturbations. Almost 80 years after the pioneering dis- 
tance - redshift diagram by Hubble, the old cosmographic 
approach appears nowadays as a precious observational 
tool to investigate the new developments of cosmology. 
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APPENDIX A: DISTANCE FORMULAE 

We derive here some useful relations for distance re- 
lated quantities as function of the redshift z and the cos- 
mographic parameters. Using their definitions in Eqs.([T]), 
it is easy to get for the fifth order Taylor expansion of 
the scale factor : 



a{to) 2 



+ ^H^t-t ) 5 + O[(t-t f 



(Al) 



with t the present day age of the universe. Note that 
Eq. (|Al[) is also the fifth order expansion of (1 + z) _1 , 
being the redshift z defined as z = a(to)/a(t) — 1. The 
physical distance travelled by a photon that is emitted 
at time and absorbed at the current epoch to is 



D = cldt = c(t - U) 

so that inserting £* = t — — into Eq. (|Aip gives us an 
expression for the redshift as function of to and D/c, i.e. 
z = z(D). Solving with respect to D up to the fifth order 
in z gives us the desired expansion for D(z) as: 



D(z) = if {V° + V]z + V\ z 2 + Vl z 3 + V\ z 4 } 



(A2) 
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with : 



with : 



D 



Vl = -(l + q Q /2) 



n x D = 1/2 , 



T) 2 — 1 J- n 4- q ° 30 

V z - 1 + qo + y - j , 



kc 2 



vi = 



- ( 1 + fflb + fflg + fflg 



1 



SQ 



2 Jo " 12 90J ° 24 



^=24 



6 + 6q Q + j - 6 



kc 2 

ff 2 n 2 
n u 



n% = 



120 



24 + 36<zo + 6gg + 8j - s 



5fcc 2 (7 + 2 9o ) 



a a H 



VI = 



l + 2q + 3q 2 + |gg + -q* 



7 



j 2 

- g<?oJo - g9 io - jo + 
1 so Iq 

In typical applications, one is not interested in the physi- 
cal distance D(z), but rather in the luminosity distance: 



Dr = 



a(t ) 



a(t - D/c) 
or the angular diameter distance 



{aoro) 



a t - D/c) 

d a = 7-t (ao^o) 

a (t ) 



(A3) 



(A4) 



24 + 48go + 18gg + 4q Jo + 12j - 2s + 24Z 



144 



3fc 2 (15 + 10<z +jo) 



144a 2 if 2 

Expressing D into Eq. (|A5[) as function of z through 
Eq. (|A2p and inserting the result into Eq. (|A3p . one ob- 
tains the desired fifth order approximation for the Hubble 
free luminosity distance dr, — Dl(z)/(c/ Hq) as function 
of the rcdshift z : 



d L {z) = V° L z + V\ z 2 + V\ ?? + V\ z 4 + V\ z r ° (A6) 
having defined : 



with ao = a(to) and 



V\ = -- (-l + «b) , 



" t0 cdt , 
sin / — — fc = 1 

'to-D/c a(t) 



ro(I>) = < 



t0 cdt 

to -D/c «(*) 



k = 



(A5) 



sinh / fc = — 1 

Jt -D/c a W 

Using Eq. (|Al[) . some cumbersome algebra finally gives: 



ro(D) 
D/a Q 
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HqD 



K 2 D 
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H n D 



HnD 
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K 



D 



HqD 

c 

HqD 



1 - qo - 3<7o + Jo + tttt 



Pi 



2 - 2cfo - I5q 2 — 15gg + 5j + lO&ojo + s 



24 



2fcc 2 (l + 3 g0 ) 
24ff 2 a 2 



m = 



-6 + 6g + 81g§ + 165gg + 105gg 
120 

lOjg - 27j - llOgojo - 105g 2 j 
120 

15g Q s + lisp + l 
120 

5fcc 2 (l + 8g + 9g§-2j ) 
120agi7 2 
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Finally, a similar procedure gives the following approx- 
imation for the Hubble free angular diameter distance 
cIa{z) = Da(z)/(c/H ) to fifth order in z : 



2fcc 2 (5 + 3g ) 



d A {z) = V° A z + V\z 2 + V\ z 3 + V\z 4 + V\ z 5 (A7) 



having set : 



2^ = 1, 



11 + 7q + 3q - j - 



kc 2 



H a 



■Da 



50 + 46 g0 + 39gg + 15gg - 13 j - lOgojo - s 
24 



_4 274 + 326g + 411g2 + 315g3 + 105 9 4 

A 120 

lOjg - 137j - 210gojo - 105g 2 j - 15go£o - 21g - / 

120 

5fcc 2 (17 + 20 g o + 9g 2 -2jq) 
120a 2 iJ 2 



Using such expressions (for k — since we have assumed 
a flat universe in the text), it is then straightforward to 
compute the quantities entering the Fisher matrix so that 
no numerical integration and differentation are needed. 
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